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ABSTRACT 

In direct sequence-code division multiple accesses systems (DS-CDMA), the signals are 
transmitted over multipath channels that introduce fading. Multipath fading along with 
multiple access interference and inter-symbol interference degrades the system performance. 
Further, simultaneous presence of multipath fading and shadowing leads to worsening of 
wireless channels. Moreover, experimental results have confirmed the presence of impulsive 
noise in wireless mobile communication channels. This paper presents a robust multiuser 
detection technique for DS-CDMA systems over generalized-K (GK) fading channels in the 
presence of impulsive noise. Robust M-decor relating detector to demodulate binary phase 
shift keying (BPSK) signals is implemented by using the particle swarm optimization (PSO) 
algorithm. The PSO algorithm is used to obtain optimum M-estimates by minimizing an 
objective function, which is a sum of less rapidly increasing function of residuals. Simulation 
results are provided to show the evolutionary behavior and efficacy of the proposed M- 
estimator in comparison with linear decor relating detector, the Huber and the Hampel 
estimator based detectors. 

Keywords: CDMA, influence function, maximal ratio combiner, M-estimator, multiuser 

detection, Nakagami fading, particle swarm optimization. 



I. INTRODUCTION 

Recent research has explored the potential benefits of 
evolutionary optimization algorithms and their 
application to multiuser detection (MUD) for direct- 
sequence code division multiple access (DS-CDMA) 
systems [l]-[4]. An adaptive robust MUD technique 
for CDMA by implementing Huber' s M-estimator 
using genetic algorithm (GA) is presented in [3] and 
shown that it is robust against heavy-tailed 
impulsive noise. Recently, particle swarm 



Optimization (PSO) algorithm has been applied for 
the MUD [4] to detect received data bit by optimizing 
an objective function. 

The generalized-K (GK) fading channel has been 
received considerable attention as it can provide a 
good fit to different fading environments such as 
Nakagami-m and Rayleigh-Lognormal [5]. 
Experimental results have confirmed the presence of 
heavy-tailed impulsive noise in outdoor mobile 
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communication channels, in radar and sonar systems 
and in indoor wireless communication channels [6] . 
Hence, this paper presents the implementation and 
performance analysis of proposed based M-estimator 
[7], [10], which performs well in the heavy-tailed 
impulsive noise, using PSO algorithm. 

II. SYSTEM MODEL 

Consider an L-user synchronous CDMA system, 
where each user transmits information by 
modulating a PN sequence over a single-path GK 
fading channel. The received signal over one symbol 
duration can be modeled as [8] 



L M-l 

r(0 = SR^Z I *l(0«/(0e 

1=1 i=0 



sfi-iTs-zM + nit) (1) 



where denotes the real part, M is the number of 

data symbols per user in the data frame of interest, T s 
is the symbol interval, «/(0 is the time-varying 

fading gain of the 1 th user's channel, (fr(t) is the time- 
varying phase of the 1 th user's channel, (z) is the i th 

bit of the 1 th user, s^(t) is the normalized signaling 
waveform of the 1 th user and n(f) is assumed as a 
zero-mean complex two-term non-Gaussian noise 
[9]. 

For synchronous case (i.e., ^ = r 2 = ... = 77 =0), 

assuming that the fading process for each user varies 
at a slower rate that the magnitude and phase can 
taken to be constant over the duration of a bit, the 
received signal can be expressed in matrix notation 
as [9] 

r(i) = AG(i) + w(i) (2) 
where r(i) □ [r\(i),...,r N (i)~] , 

T 

w(i)U[_Wi(i),...,w^f(i)] and 

1 T 

0(i)U -j^[bi(i)gi(i),...,b L (i)g L (i)] . Here, w n (i) is a 

sequence of independent and identically distributed 
(i.i.d.) complex random variables whose in-phase 
and quadrature components are independent non- 
Gaussian random variables, gi (i) is the 1 th channel 
fading coefficient and 
AU [a l ,a 2 ,...,a L ]witha l □ [aj,^,...,^]. 



It is assumed that the signal of each user arrives at the 
receiver through an independent, single-path fading 
channel. For the shadowed fading channels, a t (i) are 

i.i.d. random variables with GK distribution given by 
[12] 
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where, m is the Nakagami fading parameter that 
determines the severity of the fading, //represents 

the shadowing levels, Q 0 is the average SNR in a 

shadowed fading channel, is the modified 

Bessel function and T() is the Gamma function [9] 
In M-estimates, unknown parameters 6^ , 6 2 , ...6_ L are 
solved by minimizing a sum of function, p(-) of the 
residuals [9] 
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where p{-) represents a specific penalty function that 

is symmetric positive-definite with a unique 
minimum at zero, and is chosen to be less increasing 

than square, r n (i) and 0j (z) are the n th and 1 th 
elements of the vectors r(i) and 0{i) respectively, 

[A] n i is the nl th element of the matrix A , and 

denotes imaginary part. This paper considers 

the modified-Hampel based M-estimator to 
implement an M-decorrelating detector with penalty 
function [7] 
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where a, b and d are constants that depend on the 
robustness of the estimator [7]. 
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III. PSO BASED M-DECORRELATING DETECTOR 

PSO is a swarm intelligence method for global 
optimization modeled after the social behavior of 
bird flocking and fish schooling [4]. In PSO 
algorithm the solution search is conducted using a 
population of individual particles, where each 
particle represents a candidate solution to the 
optimization problem (3). Each particle keeps track 
of the position of its individual best solution (called 

i L \ best r best best-, j l1 11 i i i 

aspbest), =iP di ,...,p d J and the overall global 
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gbest), 



best r best best n 1 • r 1 1 .1 

g = [g J among p bests of all the 

particles in the population represented 
as p^ = [p lt d , p d ] , where p \ is the d th particle in the 

it th iteration, d = 1, 2, N P , it = 1, 2,..., Nit, N P is the 
number of particles, and Nit is the maximum number 
of iterations. Corresponding to each position, the 



it+l it it+l /r7 \ 

=Vd+ v d ( 7 ) 
where c x and c 2 are the acceleration constants 

representing the weighting of the stochastic 
acceleration terms to pull the particle to pbest and 

gbest. ju x and ju 2 are random numbers that are 

uniformly distributed between 0 and 1. Particle 
position is updated according to (6). Particle velocity 
is limited by the maximum 

i . , max r max max n 

velocity v =[v 1 ,...,v 2 J. 

Step 6: The individual best particle position is 

ifF(Pd)±F{pT) 



updated by following rule: 



„ 7 best it 

then p d = p rf 



particle velocity is v" d = [vj , v* ] [2], [4]. The steps Step ? . g ben ig ^ globa] begt position among 



involved in the PSO based M-decorrelating detector's 
implementation are [2], [3], [4]: 

Step 1: Compute the decorrelating detector output, 

Q°=RT 1 A T r. Here, r(uA t A/n) is the 

normalized cross-correlation matrix of signature 
waveforms of all users. 
Step 2: Initialization: The output of decorrelating 

detector is taken as input first particle = 0^ . 

Step 3: Fitness evaluation: The objective function (3) 
is used to find the fitness vector by substituting 

residuals. Local best position p/^is recorded 

by looking at the history of each particle and the 

particle with lowest fitness is taken as g of 

the population. 
Step 4:Update the inertia weight by using the 

decrement function w lt = (3w lt , where /?<lis 

the decrement constant. 
Step 5: Update the particle velocity by using the 
relations 



all the individual best particle positions \) d at 

the it th iteration such that F (g best ) < F (p^ ) . 

Step 8: The above steps are repeated until the 
maximum number of iterations has been 
reached. 

The computed g best value is used to evaluate average 

probability of error of BPSK demodulator, over GK 
fading channel, by an approximate expression given 
by 
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where m is the Nakagami fading parameter that 
determines the severity of the fading, ju represents 

the shadowing levels, £2 0 is the average SNR in a 



shadowed fading channel F = 



T(m)T(ju) 
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and (•) is the Whittaker function [13]. 
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TABLE I. PSO PARAMETERS USED FOR SIMULATION 



Ammeter 


V^alue 


Number of particles, N P 


20 


Maximum number of iterations, Nit 


100 


Acceleration constants c x and c 2 


2 


Maximum velocity of the particles 
v z max for all users 


2 


Initial inertial weight, w 


1 


Decrement constant, j3 


0.99 



IV. SIMULATION RESULTS 

In this section, the performance of M-decorrelating 
detector is presented by computing (7) for for 
different values of fading parameters and shadowing 
levels with least-squares, Huber, Hampel and 
modified-Hampel penalty functions. The PSO 
parameters used for simulation are presented in 

Table I. The value of g best is computed for specified 

penalty functions and is used to compute the average 
probability of error (7) of BPSK signals. In Fig. 1 and 
Fig. 2, the average probability of error versus the 
signal-to-noise ratio (SNR) corresponding to the user 
1 under perfect power control of a synchronous DS- 
CDMA system with six users (L = 6) and processing 
gain, N = 31 is plotted for moderate impulsiveness 
( s = 0.01) of noise, m = 1, \i = l. Similarly, the average 
probability of error is plotted in Fig. 3 and Fig. 4 with 
highly impulsive noise (s= 0.1) and m = 1, \i = 1. 
These simulation results show that the decorrelating 
detector with proposed estimator based detector 
performs well even in highly impulsive noise when 
compared to least-squares, Huber and Hampel 
estimator based detectors. 

V. CONCLUDING REMARKS 

Multiuser detection technique for DS-CDMA systems 
over GK fading channels using PSO was presented in 
this paper. An objective function, which is a sum of 
less rapidly increasing function of residuals, was 
used to obtain optimum estimates. An M- 
decorrelator is implemented with different influence 
functions. Simulation results shows that the 
proposed M-decorrelator performs better when 



compared to least-squares, Huber and Hampel 
estimator based detectors. 
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Fig.l Average probability of error versus SNR for user 1 for linear 
multiuser detector, minimax detector with Huber, Hampel and 
proposed M-estimator in synchronous CDMA channel with 
impulse noise, N = 31, £ = 0.01, m = 1, u = 1. 



ep>.ik<ii-0.01. k*pp«-100. L-6 

: 



- I .ca*.t Squares 

- Huhor 
-* Hanipcl 

— Proposed 




12 14 

SNR in dB 



Fig.2 Average probability of error versus SNR for user 1 for linear 
multiuser detector, minimax detector with Huber, Hampel and 
proposed M-estimator in synchronous CDMA channel with 
impulse noise, N = 31, £ = 0.01, m = 2, u = 2. 



10 



cpMnrfl I. kappa- 100. I -6 



10 



c 



10 




m- I, ft — I 



X V 
S\R in dH 



10 



12 



13 



Fig.3 Average probability of error versus SNR for user 1 for linear 
multiuser detector, minimax detector with Huber, Hampel and 
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proposed M-estimator in synchronous CDMA channel with 
impulse noise, N = 31, £ = 0.1, m = 1, u = 1. 
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Fig.4 Average probability of error versus SNR for user 1 for linear 
multiuser detector, minimax detector with Huber, Hampel and 
proposed M-estimator in synchronous CDMA channel with 
impulse noise, N = 31, £ = 0.1, m = 2, u = 2. 
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